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The fragment production cross sections for 140 MeV/nucleon **''"^Ca+^Be reactions have been 
calculated by the Heavy Ion Phase Space Exploration (HIPSE) model. Isoscaling behavior is ob- 
served. The isoscaling parameters a and /3 for both heavy and light fragments from HIPSE model 
calculations are in good agreement with the experiment data. Studies show that the potential pa- 
rameters in the HIPSE model have very small effect on the isoscaling parameters. The effect of the 
excitation energy and evaporation on the isoscaling behavior has also been discussed. 
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Motivated by the importance of the symmetry term 
in the nuclear equation of state and astrophysical ap- 
plications, interest in the isospin effect has considerably 
increased. In a series of recent papers, the scaling prop- 
erties of cross sections for fragment production with re- 
spect to the isotopic composition of the emitting sys- 
tems were investigated by M. B. Tsang et. [ll-I3|- Iso- 
topic scaling, also termed as isoscaling, has been shown 
to be a phenomenon existing in many different types of 
heavy-ion reactions [^-Q. It is observed by comparing 
the same fragments in two similar reactions that differ 
mainly in isospin asymmetry. More precisely, the ratio 
i?2i = Y2{N, Z)/Yi{N, Z) is used, where Y{N, Z) is the 
isotope yield obtained in the reactions and 2 denotes the 
more neutron-rich system. In particular, if two reactions 
have approximately the same temperature but different 
isospin asymmetry, such ratio exhibit an exponential de- 
pendence on the neutron number N and atomic number 
Z of the following form: 

R2i{N, Z) = Cexp(aiV + PZ) (1) 

where a and /3 are the isoscaling parameters and C is a 
normalization constant. In the grand canonical limit, a 
and P will have the form, a = A/i„/T and /? — Afi^/T 
where and A/itj, are the difference of the neutron 

and proton chemical potentials for two reactions. T is 
the temperature of the system in MeV. As a matter 
of fact, these equations were first obtained in the ex- 
panding emitting-source (EES) modelQ. This behav- 
ior is attributed to the difference of isospin asymmetry 
between two reaction systems in similar temperature. 
The isoscaling has been extensively studied in different 
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theoretical frameworks, ranging from dynamical models, 
such as the antisymmetrical molecular dynamics (AMD) 
model Q, to statistical models, such as the statistical 
multifragmentation (SMM) model (^]. These investiga- 
tions of isoscaling are mainly focused on light fragment. 
A few studies on the heavy projectile like residues in 
projectile fragmentation and fission fragment have been 
reported I9l-[ll|. In the studies by the statistical abrasion 
(STA) model[ll|, the isoscaling parameters decrease with 
the decrease of Z or and they are close to zero when 
Z or N become small. The decreasing trend is consistent 
with experimental data, but the value of the isoscaling 
parameter is much smaller than the data for small Z or 
N number. This is due to the simple and geometrical as- 
sumptions used in the model. For better understanding 
of the isoscaling behavior for both light and heavy frag- 
ments, the Heavy Ion Phase Space Exploration (HIPSE) 
model[l2| will be used in this paper. HIPSE parameter's 
dependence on the isoscaling parameters and extraction 
of the symmetry energy coefficient from the isoscaling 
parameters will also be investigated. 

Based on the sudden approximation and geometri- 
cal hypothesis, HIPSE model can conveniently simulate 
heavy-ion interactions at all impact parameters and thus 
can constitute a valuable tool for the understanding of 
processes such as neck emission or multifragmentation in 
peripheral and/or central collisions fl^]. Accounting for 
both dynamical and statistical aspects of nuclear colli- 
sions, HIPSE considers three different stages of the re- 
action: the approaching phase, the cluster formation, 
and the secondary deexcitation process. For the target 
and the projectile nuclei, classical two-body dynamics of 
the center of mass is assumed in the entrance channel. 
At small distance, a free parameter (aa) which defines 
the hardness of the potential has been introduced. At 
the minimal distance of approach, nucleons in each nu- 
cleus are sampled according to a realistic zero temper- 
ature Thomas-Fermi distribution. The participant and 
spectator are obtained using simple geometrical consid- 
erations. Nucleons outside the overlap region define the 
Quasi-Projectile(QP) and Quasi- Target(QT) spectators. 
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Yield ratios R2\{N,Z) for hot (lower panels) and cold (upper panels) fragments from 140 MeV/nucleon **'*°Ca+®Be 



FIG. 1 

reactions versus A'^ for the selected isotopes (left panels) and Z for the selected isotones (right panels 
isotones are shown by different symbols. The dashed and solid lines are drawn to guide the eye 
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Then two physical effects, direct nucleon-nucleon colli- 
sions and nucleon exchange, are treated in a simple way. 
The nucleons in the overlap encounter a percentage Xcoii 
of nucleon-nucleon collisions. Nucleon exchange is intro- 
duced by assuming that a fraction xtr of the nucleons 
coming initially from the target (or projectile) and be- 
longing to the overlap region are transferred to the pro- 
jectile (or target). After these preliminary steps, clus- 
ters are formed using a coalescence algorithm and propa- 
gated. To incorporate the physics of low energy reaction 
(below the Fermi energy), two nuclei can fuse. The prop- 
erties of fusion system are calculated if two fragments 
can't separate because their relative energy is lower than 
the fusion barrier. This feature leads to large Final State 
Interaction (FSI). Once all the FSIs are processed, the 
nuclei can't exchange particles anymore and a chemi- 
cal freeze-out is reached. After the chemical freeze-out, 
the total excitation energy can be determined and shared 
among fragments. At this stage, the partition is ready for 
the after-burner phase which consists of fragments' prop- 
agation in the overall coulomb field and secondary decay. 
The secondary decay is achieved using the SIMON event 
generator [l^. HIPSE model only has three adjustable 
energy dependent parameters(aa, Xtr and Xcoii)- Though 
HIPSE model is not developed to describe the projec- 
tile fragmentation process, the agreement of the fragment 
production cross sections between the calculated results 
and experimental data is quite reasonable 12 1. 



To study the isoscaling effect in projectile fragmenta- 
tion, we have chosen the reaction systems of '^^'""'Ca-f-^Be 
at 140 Mev/nucleon. Values of HIPSE's three energy de- 
pendent parameters Ua — 0.55, — 0.09, Xcoz; = 0.18 
were chosen for beam energy of 140 Mev/ nucleon [l^ . 
The production cross sections for both hot (before evap- 



oration) and cold (after evaporation) fragments are ob- 
tained directly from HIPSE model calculations. We ex- 
tract the yield ratio R2i{N, Z) from the calculated yields 
Yi{N,Z){i = 1,2) where index 2 refers to the more 
neutron-rich system. Fig. [1] shows the yield ratios of hot 
projectile-like fragments as a function of neutron num- 
ber N for selected isotopes (left lower panel) and proton 
number Z for selected isotones (right lower panel) of the 
''^''"'Ca-l-^Be reaction systems. The corresponding ratios 
for the cold projectile-like fragments are also shown in 
Fig. [T] (upper panels). In these figures, different isotopes 
and isotones are shown by alternating filled and opened 
symbols. The dashed and solid lines are just drawn to 
guide the eye. 

From the panels of Fig. [1] (hot and cold fragments), 
we observe that the ratio for each isotope Z or isotone 
N exhibits a remarkable exponential behavior. For each 
isotope (isotone), an exponential function of the form 
Cexp(aA^)(Cexp(/3Z)) was used to fit the data. By fit- 
ting the calculated points, the isoscaling parameters a 
and (3 of hot and cold fragments are obtained for all iso- 
topes {Z = 5 — 19) and isotones (A^ = 6 — 19). Using 
the same method, we have also extracted the isoscal- 
ing parameters from experimental fragment production 
cross section data[l5| and STA model calculations. The 
isoscaling parameters from experimental data and cold 
fragments calculated by HIPSE and STA model are pre- 
sented in Fig. [5] For cold fragments calculated by STA 
model, the increase of a and (3 with increase of Z and N 
is predicted. Similar trend is observed in the isoscaling 
parameters extracted from experimental data. But the 
value of the parameter from STA model is much small 
than experimental result, especially for light elements. In 
STA model, the isoscaling parameters of light fragments 
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FIG. 2: Isoscaling parameters a as a function of Z (up- 
per panel) and |/3| as a function of A'^ (lower panel) for 140 
MeV/nucleon ''^'^''Ca-l-^Be reactions. Solid symbols repre- 
sent the experimental results with the cross section data taken 
from Ref. [l^. Open squares represent calculations by HIPSE 
model. 

is almost constant for different isotopes because the exci- 
tation energy is almost same for all light fragments in the 
process of multifragmentation[TT']. In Fig. [51 a and |/3| of 
the cold fragments calculated by HIPSE model show an 
increasing trend with the increase of Z or iV. From the 
comparison of the parameters extracted from HIPSE and 
STA model with experimental results, it is demonstrated 
that HIPSE model can reproduce the isoscaling parame- 
ters of the experimental data for both projectile-like and 
light fragments quite well, especially for isoscaling pa- 
rameter \/3\. 

Isoscaling has been extensively examined in different 
theoretical frameworks, ranging from dynamical models 
to statistical models. In order to extract symmetry en- 
ergy information from isoscaling, detailed understanding 
of the sequential decay effects on the isoscaling param- 
eters is necessary. In Fig. [3] the sequential decay ef- 
fects of different model are show in one figure. There 
are four statistical models, including statistical multifrag- 
mentaton models such as SMM95, ISMM, microcanoni- 
cal multifragmentation model(MMM) and the expanding 
emitting source (EES). The dynamical models used are 
the asymmetrized molecular dynamical model (AMD), 
the Boltzmann-Nordheim-Vlasov model (BNV) and the 
isospin quantum molecular dynamical model (IQMD). 
As discussed in Ref. [16], sequential decays effects reduce 
the a values in most dynamical models and the reduc- 
tion are larger than 50%. While the effects are small in 
statistical models and its trend is not very clear. For 
140 MeV/nucleon ^*'^°Ca-f^Be reactions calculated by 
HIPSE model, the a and \/3\ from hot fragments are much 
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FIG. 3: Effect of sequential decays on the isoscaling parame- 
ter, Q!(final) and a(primary) from different models are shown 
in one panel by different symbols. The solid squares are the 
calculated results of HIPSE model, other data are taken from 
Ref. [idl and references therein. 



smaller than the parameters from cold fragments. Se- 
quential decays effects increase the parameter's values in 
the model, especially for heavy fragments. This trend 
is different from those observed in dynamical and sta- 
tistical calculations [H, [13] ■ The secondary deexcitation 
of fragments is described by the SIMON code in HIPSE 
modcl[i2]. In recent paperfli^, GEMINI IS^ replaces SI- 
MON as secondary decay code in HIPSE model. Results 
by GEMINI reproduce the cross sections of the neutron 
deficient fragments consistently better than calculations 
using SIMON. The difference of sequential decays effect 
between HIPSE and other models as shown in Fig.[3]may 
be mainly due to the evaporation program SIMON used 
in this model. A detailed study of the evaporation pro- 
cess between SIMON and GEMINI, is necessary for a 
comprehensive understanding of sequential decays effect 
on isoscaling parameters in HIPSE model. 

As a phenomenological model, HIPSE contains many 
input parameters. But for most of them, a fixed con- 
stant value is used in the model. The HIPSE model only 
has three adjustable parameters. For the reaction sys- 
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Ca-|-^Be in this paper, these parameters have 
been adjusted by comparing with experimental data and 
extrapolated with a simple function to the beam en- 
ergy of 140 Mev/ nucleon[l^ . It is necessary to study 
the HIPSE's parameter dependence of isoscaling param- 
eters. In Fig|3](from upper to bottom panel), we present 
the results by adjusting one of the HIPSE parameters 
Xtr (upper panel), XcoU (middle panel) and aa (lower 
panel) while the other two HIPSE parameter's values are 
fixed. The isoscaling parameters a show the same trend, 
and their values don't change with the adjustment, for 
most fragments. In Fig. |4l all HIPSE parameters just 
exhibit small influence on the isoscaling parameter for 
heavy fragment. So HIPSE's parameters dependence of 
isoscaling parameter can be negligible in our studies. 

In the grand-canonical approximation, assuming same 
temperature T, the isoscaling parameters a and /3 are 
proportional to the difference of the neutron and pro- 
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FIG. 4: HIPSE model parameters' dependence of isoscal- 
ing parameter a for 140 MeV/nucleon **'^''Ca+^Be reactions. 
From top to bottom, we present a as a function of Z by in- 
dependently adjusting HIPSE parameters xtr (upper panel), 
aJcoH (middle panel) and (lower panel). Experimental data 
(filled square) with the cross sections taken from Ref. [1^ are 
shown in the figures too. 

ton chemical potentials for two systems. The particular 
interest comes from their connection with the symme- 
try energy coefficient. From the difference of chemical 
potentials within the grand-canonical approximation, it 
is possible to connect the isoscaling parameter a to the 
difference of asymmetry {Zj A) between the two systems 
considered, the values of symmetry energy coefficient and 
temperature jl7| . The relation is 



a 



(|)^] = ^A[(|)2] 



(2) 



where a is the isoscaling parameter extracted from the 
calculated yields of fragments, T is the temperature, 
{Z/A)i is the ratio between nuclear charge and mass 
number. The symmetry energy coefficient Cgym is di- 
rectly related to the symmetry energy (per nucleon) of a 
given fragment having asymmetry (/ = (TV — Z)/{A)) 
In HIPSE model, the total excita- 
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tion energy is determined by the total energy balance 
and shared among fragments. In upper panel of Fig. [Sj 
we present the extracted parameter a as a function of 
excitation energy {E* /A), a show a decreasing trend 
with the increasing of E* /A. According to Eq. ([2]), if we 
use the Fermi-gas relationship between excitation energy 
per nucleon and temperature E* /A = ^T^ to calculate 



T, with the inverse level density parameter a = 8 — 13 
(a = 10 is used in this paper), the symmetry energy co- 
efficient {Csym) could be extracted. Results extracted 
from a are shown in the lower panel of Fig. [SJ The sym- 
metry energy coefficient from a is around 20 Mev and 
decreases with the increasing of excitation energy E* /A. 
Similar dependence was also observed in the experimen- 
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FIG. 5: The isoscaling parameter a (upper panel) and the 
extracted symmetry energy coefficient from a {Csym, lower 
panel) as a function of excitation energy (E* /A) for 140 
MeV/nucleon ''*'''°Ca-l-^Be reactions. The solid squares are 
the experimental data taken from Ref. 0]. 

tal dataQ but their Csym values is a little bit smaller 
than HIPSE's result. The experimental studies of heavy- 
residue isoscaling reveals the gradual decrease of the sym- 
metry energy coefficient with increasing excitation energy 
in the rang eE*/A = 2.0 - 2.9 0,[il, a similar trend is 
obtained by HIPSE in a big range E* /A = 2.0 - 4.5. 

In conclusion, the fragment production cross sections 
of 140 MeV/nucleon *®'*'^Ca-|-^Be reaction systems have 
been studied by HIPSE model. Isoscaling behavior is 
observed in the ratio of the fragment yields from the 
two systems. The extracted isoscaling parameters from 
HIPSE model calculations are in good agreement with re- 
cent experimental data for both projectile-like and light 
fragments. The excitation energy dependence of sym- 
metry energy coefficients obtained from the calculations 
is also consistent with the experimental data. Sequen- 
tial decays and parameters of HIPSE model have some 
effects on isoscaling parameters, especially for heavy frag- 
ments. Thus HIPSE model could be used to investigate 
the isoscaling phenomena of fragments produced in cen- 
tral and peripheral collisions. 
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